ID Quantum Liquids with Power-Law Interactions: the Luttinger Staircase 
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We study one dimensional fermionic and bosonic gases with repulsive power-law interactions 
with /3 > 1, in the framework of Tomonaga-Luttinger liquid (TLL) theory. We obtain 
an accurate analytical expression linking the TLL parameter to the microscopic Hamiltonian, for 
arbitrary /3 and strength of the interactions. In the presence of a small periodic potential, power- 
law interactions make the TLL unstable towards the formation of a cascade of lattice solids with 
fractional filling, a "Luttinger staircase" . Several of these quantum phases and phase transitions are 
realized with groundstate polar molecules and weakly-bound magnetic Feshbach molecules. 
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There is presently considerable interest in quantum de- 
generate gases with long range interactions in reduced 
geometries [ij. This is motivated by recent experiments 
with polar molecules 0], where electric dipolc moments 
associated with rotational excitations lead to strong, 
anisotropic dipolar interactions but also by experi- 
ments with atomic gases with strong magnetic dipoles Q ■ 
For polar molecules, electric dipoles can be manipulated 
with external microwave AC and DC electric fields, which 
provides a toolbox to tailor the many-body interactions, 
and in combination with optical trapping in ID or 2D 
promises the realization of stable exotic, strongly corre- 
lated quantum phases with long range interactions . 

An intriguing example is given by polar molecules 
trapped in a ID wire [see Fig. [U^a)] [gHsI, where long 
range interactions compete with an optical lattice in a 
commensurate - incommensurate transition. In the zero- 
tunneling limit in a deep lattice, this leads to the forma- 
tion of a devil's staircase, that is, a continuous and non 
differentiable (Cantor) function for the ground state fill- 
ing fraction as function of the chemical potential ^, stud- 
ied first in the context of atomic monolayers adsorbed on 
solids While recent studies at finite hopping have 

already shown a modification of this structure in a deep 
lattice [m , the challenge is now to investigate the quan- 
tum regime where large kinetic energies compete with 
both interaction strengths and periodic confinement. 

Below we show that using bosonization techniques 
the classification of quantum phases can be derived ana- 
lytically for all power law interactions Ci3/\x\^ with /3 > 1 
and for arbitrary relative strengths of the kinetic energy 
and the long-range repulsion. Remarkably, the parame- 
ters of the effective bosonized theory can be accurately 
obtained in analytical form for all /? in terms of the mi- 
croscopic Hamiltonian, even in the absence of integra- 
bility. This provides us with a universal phase diagram 
where the cases of repulsive Van der Waals (/? = 6) and 
dipolar interactions (/5 = 3) should be accessible in po- 
lar molecule experiments [5|. In contrast to the classical 



devil's staircase, where lattice solids are stable over a 
finite interval in /i for every rational filling fraction be- 
tween and 1, and the total measure of such interval 
exhausts the full range of ji, we find that in the ID quan- 
tum case large kinetic energies prevent the formation of 
ordered states, where the average interparticlc distance 
is not constant. This drastically reduces the number of 
"steps" in the staircase to a number not dense in the 
interval ]0, 1], i.e. a Luttinger staircase. Signatures of 
these quantum phases are excitations in the form of soli- 
tons and breathers, detectable via Bragg scattering. 

We assume that the polar molecules are polarized by 
external electric fields, and confined to a ID geometry, 
e.g., by a sufficiently deep 2D optical lattice with fre- 
quency Lo^. The shape of the long distance interactions 
can be tuned by coupling the lowest rotational mani- 
folds of each molecule with DC and microwave AC fields. 
As shown in Refs. [sj we can tune between /? = 3 with 
C3 = (P /to, where d is the dipole moment induced by an 
electric field i^DC and eo the vacuum permittivity, and 
/3 = 6 with Cq oc d'^/hA, where A the detuning of a mi- 
crowave field Eac coupling the ground to the first excited 
rotational manifold. For average interparticlc distances 
a ^ {Cp/hwi_Y/^ the gas dynamics is one-dimensional 
and microscopically described by the Hamiltonian 
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Here, ipi^) is a field operator for molecules, which can 
be either fermionic or bosonic, m is the particle mass, 
and U{x) = U sin^(27ra;/A) is a weak periodic poten- 
tial, as provided by a shallow optical lattice of strength 
Ul = U/Er < 1, with Er = h^/2m\^ and A the lattice 
wavelength. For the case /? = 3 with experimentally rel- 
evant molecules such as LiCs, RbCs or KRb molecules 
(rfmax = 5.6, 1.25 and 0.5 Debye, respectively) and con- 
finement Lo±_ = 27r X 100 kHz, {Cp /ftwj^)^/^ is of the order 
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of 360, 130 and 80nm, respectively (HI. [l^. 

In the absence of an optical lattice (Ul = 0) the short 
range character of power-law interactions with (3 > 1 
allows a description of the low energy physics in terms of 
Tomonaga-Luttinger liquid (TLL) theory [l3 , 14 1 . Here, 
we first consider the bosonic case, and then discuss the 
differences with the fermionic one. The TLL effective 
Hamiltonian is given by 0, EE] 

H = hv fdx [id^9{x))yK + Kid^(f>{x))^] /(27r). (2) 



Here, the field i^{x) in Eq. ([T]) is replaced by ip{x) ~ 
\/n + dxd{x) /tt exp[i(j){x)] in a hydrodynamic approach, 
and dxO{x) and dx4>{x) characterize the long-wavelength 
fluctuations of the density n and of the phase (t>{x), re- 
spectively, with [dxO{x), (j){y)] ~ in6(x — y). The liquid 
is completely characterized by the sound velocity v and 
Luttingcr parameter K = hmT/{mv), which determines 
the algebraic decay of the correlation functions 



{n{x)n{x')) ^ \x — x' 



-2^, {^{x)i}^ {x')) ~ \x-x'\-^'^'^. 

(3) 

In general, K can be related to the microscopic param- 
eters of the Hamiltonian only for exactly solvable models, 
e.g., contact interactions or /3 = 2 [Calogero-Sutherland 
(CS) model]. Below we show that the dependence of K 
on the microscopic parameters in ([T]) can be given ana- 
lytically for arbitrary shape and strength of interactions. 



-1/2 



(4) 



with n^^'^Rp the dimensionless interaction strength, and 
Rp = mC p / {2'Kti?') (see Figl3]). In contrast to famil- 
iar bosonic gases with contact interactions where K > 
1 13, iBj[i3] J long-range power-law interactions constrain 
K to values 1 > A' > 0, where K = 1 corresponds to 
the Tonks- Girardeau gas limit and AT = to a system 
with long-range order [l3|. Eq. (|4]) allows us to readily 
determine the phase diagram for Ul = 0, by compar- 
ing the relative decay of the correlation functions in ([3]): 
a crossover from superfluid (SF) to charge-density wave 
(CDW) behavior takes place at Kc = 0.5. The fermionic 
gas is also described by Eqs. ^ and (|4|), however its 




FIG. 1: (color online)(a) Experimental setup (sketch): an 
array of ID polar molecular gases is formed along x (green 
tubes); molecules are polarized perpendicular to x. (bl) 
Groundstate configuration in the solid phase with filling 1/p — 
1/3. (b2) soliton and antisoliton excitations with repulsive in- 
teractions with 1/p = 1. (b3) a breather, 1/p = 1 (text). 
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FIG. 2: (color online) (a) Commensurate phase diagram for 
bosons with dipolar interactions /3 = 3, and lattice depth 
Ul ~ 0.1. Physical configurations correspond to commen- 
surate fillings nX/2 = 1/p, with p G N (horizontal lines are 
guides to the eye for p < 10). Quantum phase transitions 
from a TLL to a lattice solid [or Mott insulator, (MI)] occur 
for each 1/p at the position of the dots on the continuous line, 
while red and blue dots on dashed lines indicate crossovers. 
Mil and MI2 indicate MI with solitonic and breather excita- 
tions, respectively (see Fig.l). (b) Phase diagram at commen- 
surate filling 1/p = 1/3 in the Ul vs uR^ plane. Continuous 
line: quantum phase transition between a TLL and a lattice 
solid. The phase diagram for fermions is identical to the one 
for bosons, except the TLL is always a CDW. 



phase diagram displays a CDW behavior at all interac- 
tion strengths. In addition, correlation functions in ([5]) 
have a slightly different long-distance decay [l^. In the 
following, statistics will not be relevant, and thus we deal 
with both cases at the same time. 

The Luttingcr staircase: A cascade of insulating lat- 
tice solids can be realized from a TLL with power-law 
interactions, by introducing a vanishingly-small periodic 
lattice potential, as provided by the shallow optical lat- 
tice U{x) in Eq. ([T]). Combining the complete density 
operator with the periodic term U{x) [l7| . one obtains 



^ [/^ = ^ I j cos[r2e{x) + Qrx\d2 



(5) 



with Qr{n,X) = 27r (rn - 2/A) and r e N. Here, U = 
7r7iC/A^/2, where A is a cutoff that fixes the energy scale 
of the initial Hamiltonian [l^. The term Qr in Eq. ([5|) 
is responsible for a possible competition between two 
lengths scales: the interparticle distance 1/n and the lat- 
tice period A/2. We can then distinguish two different 
situations: a commensurate one with Qr = 0, where the 
length-scales do not compete, corresponding to the con- 
dition 2/ (An) = p G N, and an incommensurate one with 
Qr 7^ 0, where a competition is present. 

In all commensurate cases 2/(An) = p G N, the most 
relevant term due to the optical lattice in ([5]) in the 
renormalization- group sense is Up. Keeping only this 
term, the system becomes equivalent to a sine-Gordon 
model [l3| , where the scaling dimension of the cosine op- 
erator is affected by the interparticle interactions through 
K and by the optical lattice through U. For p > 1 we 
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then expect that for weak interactions and small depth 
of the optical lattice the cosine term in ([5]) is irrelevant 
and the TLL liquid is preserved, with correlation func- 
tions decaying polynomially as in ([3]). However, when the 
cosine term is relevant, we have a non-zero expectation 
value (cos[p-\/47r0(a;)]) ^ and the system is pinned on 
the lattice. This pinning corresponds to the breakdown 
of TLL and to the formation of a lattice solid, or Mott 
insulator (MI), with particles localized at individual sites 
of the lattice, every p lattice sites. In this phase, the ex- 
citation spectrum is gapped and the off-diagonal corre- 
lation function decays exponentially. Using Berezinskii- 
Kosterlitz-Thouless (BKT) scaling near criticality [l3| . 
we find that for each p the gapped phase occurs for 



2 + Ul> Kp' 



(6) 



with Ul = U/Er < 1. Eq. ^ is remarkable in that 
it shows that power-law interactions make possible the 
realization of an infinite series of gapped phases at lattice 
filling less than one [Fig.[Hbl) for p = 3]. The case p = 1 
is peculiar as Eq. ([6|) is always satisfied, implying a lattice 
solid for a vanishingly small Ul- The cascade of solids 
with p € N corresponds to a quantum version of a Devil's 
staircase structure, where large kinetic energies of order 
of Eji prevent the formation of ordered states where the 
average intcrparticlc distance is not constant. This is in 
contrast to the classical Devil's staircase of the Frcnkel- 
Kontorova model [lo|, where commensurability is also 
allowed for rational fillings r S Q 7^ N. Evidence of this 
classical case have been recently found in the deep lattice 
limit of Refs. 
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for /? = 3, in ID and 2D. 
Equation ([6]) shows that a gap is favored by high densi- 
ties, strong interactions and finite (small) lattice depths. 



For C/l < 1 a good estimate for the gap A is given by [1 



ttK 



A = 



L{4-2K)J 



16 r[i 




(7) 



When K is close to 1/p^, A approaches the massive 
fermion limit A ~ [//2 recently observed for p = 1 with 
contact interactions [l3, , whereas close to the BKT 
transition it closes exponentially. In the vicinity of the 
BKT transition excitations are of the soliton/antisoliton 
type, which in the massive fermion limit correspond 
to weakly repulsive particles and holes. Fig. [Hb2). In 
contrast to contact interactions, power-law interactions 
allow one to tune the sign of the soliton-antisoliton 
interactions from repulsive {K > 1/p^) to attractive 
[K < 1/p^), giving rise to soliton-antisoliton bound 
states called breathers [Fig. [TJb3)]. These excitations 
are confined in space but oscillatory in time, and are 
stable solutions of the equations of motion for the sine- 
Gordon model [isj . The number of breather excitations 
is A/" = 2(1/ A' - 1), with energy 

Mbin') = 2A sin[7rn7(4/X - 2)], n' < N. (8) 
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TABLE I: Minimal lattice depth Ul(1/p) for a lat- 
tice solid at filling 1/p, for groundstate polar molecules 
(KRb, RbCs, LiCs) and magnetic Feshbach molecules 
(^^Cr2,^^*Dy2,^®^Er2); an arbitrarily small periodic potential 
pins the TLL for 0+. Lattice depths with Ul{1/p) ^ 2, where 
the sine-Gordon picture breaks down [1^ , are not considered. 



For K < l/(2p^) breathers are the lowest-energy exci- 
tations, qualitatively changing the spectrum of the insu- 
lating phase with respect to the familiar case of contact 
interactions. Strong power-law interactions will allow for 
an unambiguous observation of this localized topologi- 
cal excitations, with applications ranging from Josephson 
junctions to conjugated polymers, see below [20j . In the 
vicinity of commensurate fillings (p g N), the system ex- 
hibits a gap as long as the energy shift due to Qr 7^ 
remains small with respect to A' A/2, the energy required 
to add a particle: this ensures stability of the phases 
above with respect to small density changes. Above a 
critical Qc a commensurate-incommensurate phase tran- 
sition takes place from an insulator to a gapless phase, 
similar to contact interactions [l3|- For generic values of 
Qr^O, the TLL is stable. 

Figure [51^a) shows the commensurate phase diagram for 
the case of bosonic particles with dipole-dipole interac- 
tions /3 = 3 as a function of the lattice filling 1/p = nX/2 
and the strength of interactions R3. The lattice depth 
is Ul = 0.1. For each p E N, the BKT quantum phase 
transition occurs at the position of the dot along the 
continuous line, while dots on the dashed lines charac- 
terize crossovers. The regions denoted as Mil and MI2 
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FIG. 3: (color online) (a) TLL parameter K vs the dimen- 
sionless interaction strength n^'^Rp for dipolar interactions 
/? = 3. Line: analytic result Eq. ((4|. Squares and dots: 
quantum Monte-Carlo (QMC) results of Refs. Q and re- 
spectively. Inset: K vs n^~^Rp with 13 — 2. Continuous line: 
Eq. Q. Dashed line: exact Calogero-Sutherland model, (b) 
TLL velocity v vs nRz for (3 — 3, with no = h/(y/2mRz). Line 
and dots: analytic and QMC results of Ref. respectively. 
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correspond to MI with soliton/antisoliton and breather 
excitations, respectively, and the dashed hne signals the 
crossover for TiT = 1/ (2p^) . Panel (b) shows the transition 
between the TLL and solid behavior as a function of the 
lattice depth and nR^ for the case of p = 3. The phase di- 
agram for fermionic particles is identical to Fig. [51 except 
that the TLL phase is always a CDW. Phase diagrams 
for /3 ^ 3 look qualitatively similar to Fig. [21 

In Tabic U we list the estimated minimal lattice depth 
necessary to realize a MI with filling 1/p, for a few 
groundstate polar molecules. In addition, we report esti- 
mates for magnetic Feshbach molecules, where the ruag- 
netic dipole moment is taken as twice the atomic one j22l | . 
The realization of insulating states with, e.g., p = 2 will 
help stabilize highly-excited Feshbach molecules against 
three-body recombination, opening the way to the real- 
ization of strongly-correlated lattice phases. 

Analytical expression for K: For Ul = 0, Eq. (|T]) 
describes an effective, strictly one-dimensional, scale 



invariant theory, dependent only on 



Rfs at all 



length/energy scales, 
teractions read V{y = 



After rescaling, dimensionless in- 
rn) = Rj3n^~'^ /y^ , and ultraviolet 
divergences in Eq. (|T]) can be treated by introducing a di- 
mensionless cut-off A We choose A such that V{y) cx 
{y + A)~^, so that scale- invariance is preserved explicitly: 
the shape of the regularized potential is independent of 
and n. We can now fix A self-consistently in the 
effective ID theory, as follows [consistency with the mi- 
croscopic derivation of (|T]) implies A/n > (C^g/ftwj^)^/^]. 
We compute K analytically in the strong- {n^~^Ri3 ^ 
1) and in the weak-coupling (n^~'^Rp <C 1) limits as 
lU - 7r/[/3(/3 + l)C(/3)i?«n^-V2]^/' i and = 
(l + n'3-2i?^V^(0)/27r)-i/2 [13, respectively, with V{Q) = 
A^~^ /{P — 1) the Fourier transform of V{y) at zero- 
momentum, and C, the Riemann Zeta-function. Due 
to the similar functional dependence, we then fix A = 
[/3(/3- l)(/3-f l)C(/3)/7r]i/(i-'') by matching = for 
n^~'^Rl3 > 1 and obtain the approximate Eq. ^ |23l] . 

Expression ([4]) compares favorably with known exact 
results. In Fig. [Slja) (Inset) we compare it to the exact 
expression Kcs = 2/(1 ^/T+TR^) for ^ = 2, which 
we derive from the Bethe-Ansatz solution of the CS- 



model [21|. We find quantitative agreement between the 
two curves for the entire range of interaction strengths 
< i?2 < 100; with a maximal relative difference of 
about 5% at i?2 — 1; and recover the n-indepcndence of 
the CS model |2l| . Furthermore, in the main figure wc 
compare K = + 0.73ni?3, as derived from Eq. (|4]) 

for /? = 3, to the numerical quantum Monte-Carlo results 
of Refs. and Q (black squares and red dots, respec- 
tively), finding good agreement for < uR^ < 1000. In 
panel (b), we also plot the velocity v in the same range 
of ni?3 values, finding excellent agreement with the re- 
sults of 0- This fixes the phenomenological parameters 
in the effective Hamiltonian We are not aware of 
exact results for /? > 3 to compare with our predictions. 



The extension of the techniques described here to 
several species will enable a microscopic treatment of 
strongly correlated phenomena in mixtures of polar 
molecules in single- and multi-tube configurations, as rel- 
evant to experiments Q , in particular exotic phases such 
as bond-ordered density waves and trimer liquids. 
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